Accepted for publication in ApJ Letters on 19 December 2005 

Preprint typeset using I^T^X style cmulatcapj v. 6/22/04 



FAR-ULTRAVIOLET OBSERVATIONS OF THE NORTH ECLIPTIC POLE WITH SPEAR 

Eric J. Korpela 1 , Jerry Edelstein 1 , Julia Kregenow 1 , Kaori Nishikida 1 , Kyoung-Wook Min 2 , Dae-Hee Lee 2,3 , 
Kwangsun Ryu 2 , Wonyong Han 3 , Uk-Won Nam 3 , Jang-Hyun Park 3 

Accepted for publication in ApJ Letters on 19 December 2005 

ABSTRACT 

We present SPEAR/FIMS far-ultraviolet observations near the North Ecliptic Pole. This area, at 
b ~ 30° and with intermediate Hi column, seems to be a fairly typical line of sight that is representative 
of general processes in the diffuse ISM. We detect a surprising number of emission lines of many 
elements at various ionization states representing gas phases from the warm neutral medium (WNM) 
to the hot ionized medium (HIM) . We also detect fluorescence bands of H2 , which may be due to the 
ubiquitous diffuse H2 previously observed in absorption. 
Subject headings: ISM: general, ISM: lines and bands, ultraviolet: ISM 



1. INTRODUCTION 

The North Ecliptic Pole (NEP) is a region of the sky 
with no obviously unusual features. Many prior space- 
missions (Einstein, ROSAT, COBE, IRAS, etc.) have 
conducted deep surveys of this region. Several ground 
based surveys have also concentrated on t his region as 
a complement to the s pace-based surveys l)Labov et al.l 
Il989t lElvis et al]ll994lL Its location at moderate galac- 
tic latitude (b=+29°) places it clear of the bright stars 
of the galactic plane. In this region, the Galactic 
N(HI) varies from about 2xl0 20 cm" 2 to 8xl0 20 cm" 2 
which corresponds to an dust opacity of Td us t =0.6- 
2.3 at 1032A()Sasseen et alJ 12002ft or T dust =0.3-1.4 at 
1550A llSavage fc Mathi3ll97flh . 

The Spectroscopy of Plasma Evolution from Astro- 
physical Radiation (SPEAR) instrument, designed for 
observing emission lines from the diffuse ISM was 
launched in late 2003. SPEAR is a dual-channel FUV 
imaging spectrograph (short-A (S) channel 900 - 1150 
A, long-A (L) channel 1350 - 1750 A) with A/AA -550, 
with a large field of view (S: 4.0°x 4.6', L: 7 . 5° x 4.3') im- 
aged at 10' resolution. See lEdelstein et aD l)2005albl ) for 
an overview of the instrument, mission, and data anal- 
yses. SPEAR sky-survey observations consist of sweeps 
at constant ecliptic longitude from the NEP to the south 
ecliptic pole through the anti-solar point. The duration 
of each sweep varies between 900 and 1500 s. Each sur- 
vey scan includes the region near the north and south 
ecliptic poles resulting in large exposure times in these 
regions. We report on deep SPEAR observations of FUV 
emission spectra from a 15° radius region centered on the 
NEP. 

2. OBSERVATIONS AND DATA REDUCTION 

The NEP data analyzed here include sky-survey and 
pointed observations that occurred between 8 November 
2003 to 10 November 2004. To remove times of high 
background and bright stars from the data set, we have 
excluded times during which the count rate exceeds 100 
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s" 1 . (In comparison, a typical photon count rate is <20 
s .) Because stars transit the instrument slit in — 5 
s during survey sweeps, this results in minimal loss of 
observing time. The area of sky lost due to a star is 
similar to the instrument imagi ng resolution acro ss the 
slit width (~ 10' x 5'). See lEdelstein et all l|2005aD for a 
description of the data reduction pipeline. 

To further remove stellar contamination from our 
dataset, we excluded locally intense spatial bins (~30% 
of the viewed area), defined to have >3x the median 
count rate. This level corresponds with a slope change 
in the histogram of the log count rate, indicative of a 
difference in distribution of scattered vs direct starlight. 
Inclusion of unresolved starlight should not be a factor 
other than to raise emission-line determination errors. 
To obtain the highest sensitivity possible, we have binned 
the entire dataset into a single spectrum including a to- 
tal of 3.5 xlO 5 counts over an effective full-slit exposure 
of 13.5 ks. 

3. SPECTRAL MODELING 

The raw SPEAR spectrum is composed of many com- 
ponents: detector background, dust-scattered stellar con- 
tinua, direct and instrument ally scattered airglow, and IS 
atomic and molecular emission. To model a spectrum we 
simultaneously fit all of the components to the observed 
spectrum. 

The spectral shape is distorted by the detector elec- 
tronics and by the data processing pipeline. These 
distortions appear as fluctuations in a flat-field image. 
Any broad-band component used in our model spec- 
trum must be multiplied by this non linearity function. 
<|KorT)elaet alJl2003t iRhee et alJl200l . 

Both S and L are affected by instrumental background 
due to cosmic rays, radioactive decay within the detector, 
and thermal charged particles entering the instrument. 
These sources of background are relatively uniform across 
the face of the detectors. Although the charged particle 
rate can change with time and orbital position, the po- 
sitional distribution of these background events is domi- 
nated by static fields within the instrument and is rela- 
tively constant. A sum of many shutter-closed dark ex- 
posures multiplied by a rate factor can be used to model 
this contribution to the spectrum. 

In L, the largest broad-band spectral component is 
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Fig. 1. — S channel H2 and line emission models. Upper panel: best fit H2 fluorescence model is shown against the residuals following 
subtraction of the continuum/background model. The H2 features' positions are well predicted, but their ratios are not. The H2 is 
oversubtracted near 1100 A and undersubtracted shortward of Ly {3 A 1026. Features at A < 1000 A are underpredicted. Lower Panel: The 
best fit atomic line spectrum is plotted against the residuals following subtraction of the H2 spectrum. Probable H2 features and several 
as yet unidentified features remain throughout the band. 



dust scattered stellar continua. We model the input spec- 
trum to the dust scattering process as a sum of spectra 
for upper main sequence (UMS) stars. Each spectrum 
is weighted by using a power law UMS luminosity func- 
tion (dN/dM v oc 10 qM ", where M v is the absolute mag- 
nitude and dN/dM v is the volume density of stars per 
unit magnitude) and a dilution factor dependent upon 
the weighted mean distance of the stars. The spectral 
sum is multiplied by a dust opacity s caled to a va riable 
N(HI) and by a dust albedo function l)Drainel 2003) with 
a variable total dust surface area. The best fits for this 
function typically occur very near to the meas ured UMS 
luminosity function with a ~ 0.2 llReedl ::2001h Because 
the stellar features in this spectrum cannot be directly 
measured, there may be some bias in our measurements 
of the overlying resonance lines (Civ AA 1548,1551, Ovi 
AA 1032,1038, SiivAA 1394,1403, etc.) of order - 1000 
ph s _1 cm~ 2 sr _1 (LU). We include this effect in our 
error analysis. Our best fit continuum, which could in- 
clude unresolved stars, is 300±5q° ph s _1 cm" 2 sr _1 - 

A (CU) in L. In S, due to its lower sensitivity and 
higher non-astrophysical background, we are only able to 
place a statistical upper limit of 1200 CU. However, given 
the wavelength dependence of the dust albedo, we expect 
that the true S astrophysical continuum is no higher than 
in L. 

We first attempted to fit atomic line emission by using 
linear combinations of collisional ionization equilibrium 
(CIE) plasma models with a limited abundance gradi- 
ent vs T. We attribute this method's lack of success to 
photoexcitation of resonance lines, non-CIE effects and 
abundance variations that are not T related. To fit the 
line emission we have, instead, fit the spectrum of each 
species separately with a per-species intensity and T pa- 
rameter. The emission line spectra are convolved with 
a A-dependent Gaussian to model the instrument spec- 
tral res olution. We use th e atomic parameters of CHI- 
ANTI l)Young et al.l l2003|) to determine the spectrum 
of the species vs T. This spectrum is then normalized 
and multiplied by the intensity parameter. When CHI- 
ANTI parameters are not available for a species, we use 
a spectrum calculated using the NIST Atomic Spectra 
Database. 

We then consider the effects of self-absorption of spec- 
tral resonance lines. In the case of highly ionized species 
we assume all lines are optically thin. We also as- 
sume this to be the case for species that would be com- 



pletely photoionized in the WNM (ionization potential 
(IP) <13.6 eV). For species which are likely to be abun- 
dant in the WNM and warm ionized medium (WIM), (IP 
13.6 eV) we assume that resonance lines are optically 
thick (r >10) and that the distances to the illuminating 
stars are large compared to the r = 1 surface. In cases 
where the ground state is split, absorbed resonance line 
photons will be converted to the excited state, greatly 
changing the relative contributions in the multiplet. We 
model this by attenuating the ground state transition by 
a large factor. We will develop a more accurate way of 
treating these optical depth effects in future work. 

4. DISCUSSION 

The spectrum of the bright Lyman-series geocoronal 
emission, easily visible in Fig. ^ can be modeled using a 
temperature and total line flux. The largest S broadband 
component is scattered geocoronal Ly a A 1216 emission. 
This line occurs in the bandpass gap between S and L. 
L is immune to Ly a contamination because it includes 
a CaF2 filter which is very nearly opaque to Lya radi- 
ation. We fit the S scattering with two components, an 
exponential scatter oc exp(— |A — Xq \ /a) with a ~ 200A, 
and an isotropic component. The scattering, summed 
over the entire spectrum, is ~ 0.08% of a ^5 kRy Lya 
incident flux. Because the Lya transition is optically 
thick, and because the scattering fractions are indepen- 
dent parameters, we fix the estimated incident Lya flux 
at HOOx the Ly/3 flux which helps the solution to con- 
verge. While the fractional scatter and exponential width 
of these components are free parameters in our spectral 
model, the values we obtain vary little from one dataset 
to the next. 

Our model includes spectral features of atmospheric 
Ol, Cl, and Ni that are much fainter than when seen in 
day-time observations (e.g. FUSE), as SPEAR only ob- 
serves during orbital eclipse. We detect a spectral feature 
near 990 A that is consistent with the Ol A 989 transi- 
tion. However, this feature may be blended with Nm A 
990. If the entirety of the 990 A feature is due to Ol, 
the lack of detectable Ol A 1356 emission would indicate 
T > 10 3 ' 6 K, which is plausible for interstellar Ol. We do 
detect Ol A 1356 in data that is not count rate filtered. 
The other possibility is that the 990 A is due to Ol at 
T ~ few hundred K which seems unlikely even for atmo- 
spheric Ol. The most likely explanation for the lack of 
any detection of Ol emission at A ^ 990A is that a large 
portion, if not all, of the 990 A emission is due to Nm. 
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We detect relatively strong emission from H2 fluores- 
cence bands even in this region of intermediate N(HI). 
This is not entirely surprising since absorption due to 
H2 is a ubiquitous feat ure of FUV spectra taken by 
FUSE. (|Shull et alJl200(l We fit these features by using 
interpolation between of a sparse grid of H2 models (B. 
Draine 2004, private communication). The upper panels 
of Figs. HI show the best fit H2 model plotted against 
the residuals of the spectrum following background, con- 
tinuum, and scattering subtraction. We detect H2 emis- 
sion features in both spectral bands. In S we detect the 
H 2 AA 965 band (to the left of the HI A 972 line) and 
the H 2 AA 986 band (to the left of the Oi,Niii AA 989 
feature). The H2 A A 1053 and A A 1097 bands are not 
detected at any significance. In L, H2 provides contribu- 
tions throughout the spectrum. We attribute the triple 
peaked band structure between 1350 and 1380 A, the 
peak near 1440 A, the broad features AA 1520-1530, AA 
1570-1590, and the peak near 1608A to H 2 . Note that 
these features are not well fit by our H2 model. It is 
not surprising that our sparse linear model with a sin- 
gle dust absorption column does not fit these data since 
the observed emission derives from many lines of sight 
with varying illumination, gas temperatures and absorb- 
ing columns. In Figs. El lower panels we show the 
residual spectrum following subtraction of the H2 com- 
ponent. 

We see intriguing indications of N2 fluorescent emis- 
sion , most notably the features at 1434A, 1445A, 1475A, 
and 1660A. FUSE has detected absorption due to N2 
along some sight- lin es, albeit at shorter wavelengths 
ijKnauth et alJl200^) . We cannot yet discount the pos- 
sibility that these N 2 features are atmospheric, as near 
the magnetic poles, Nj" can be transported to high alti- 
tude. Further study to determine whether the N2 emis- 
sion varies with spacecraft orbital position or with galac- 
tic latitude of the observation should resolve this ques- 
tion. If confirmed as astrophysical, this would be the 
first detection of fluorescent N 2 emission from the diffuse 
ISM. 

SPEAR has discovered some surprising spectral fea- 
tures, such as the resonance lines of the singly ionized 
species, Sill and Alii. The Alii feature at 1671 A very 
nearly overlies the Om] AA 1665 which we do not de- 
tect in this dataset. In a lower resolution instrument 
these spectral features could be blended resulting in con- 
fusion of the Om] and Alii emission. We think it is 
likely that at least some of the prior claimed detections 
lIMartin fc Bowverlll99(il) of Om] were in fact misidenti- 
fied detections of Alii AA 1671. We attribute the feature 



TABLE 1 

Observed warm/hot ISM emission line intensities 



Species 


A 


h 


Statistical 


Model 




A 




Error 


Uncertainty 




LU a 


LU 


LU 


Om] 


1667 


<500 6 






Orv] 


1400 


1980 


±220 


±300 


Ovi 


1032,1038 


5724 


±570 


±1100 


Cm 


977 


1260 


±270 




Civ 


1548,1551 


5820 


±280 


+0 

-1490 


Sin* 


1532 


2430 


±200 


+630 
-0 


Siiv 


1394,1403 


1430 


±120 


±160 


Nil 


1085 


8500 


±730 




Nm 


990 


<2300 






Niv 


1485 


1984 


±170 


±450 


Alii 


1671 


5608 


±525 




Hen 


1640 


1850 


±180 


+0 
-500 


Nevi 


999,1005 


<4100 







Note. — a LU=ph s 1 cm 2 sr 1 b Upper limits are 90% 
confidence 



at 1533 A to the Sin* A 1533.4 transition into the up- 
per fine structure level of the ground state. The offset 
between the model feature and the observed line may be 
due to an unc orrected continuum feat ure underlying this 
line (Fig. 3 of lEdelstein et all <|2005a|) ). Since Sin is the 
predominant state of silicon in both the WNM and WIM, 
the true Sin A 1526.7 ground state transition is likely to 
be optically thick and therefore undetectable. Our spec- 
tral models are improved by the addition of Fen transi- 
tions throughout band, most notably the Fen AA 1564, 
which fills the gap between Civ AA 1548,1551 and the AA 
1570-1590 H2 band. It is also possible that there is a con- 
tribution to the H 2 AA 1600-1620 band due to the Fell 
AA 1608 resonance lines. These lines have a fairly low 
oscillator strength and are expected to be optically thin 
along this line of sight. However, given the low oscillator 
strength, we expect these lines to be faint compared to 
the Sin and Aln features. 

Emission lines of highly ionized gas expected to be gen- 
erated in a CIE plasma at HIM temperatures are shown 
in Tabled Of these we detect (> 3a) Cm A 977, Civ A A 
1548,1551, Siiv AA 1394,1403 and Ovi AA 1032,1038. The 
Civ AA 1549 and Siiv AA 1400 emission overlie prominent 
stellar absorption features which are likely to be present 
in the dust-scattered stellar radiation. Since the depth 
of these features is unknown we set the lower range of 
our model uncertainty to be the intensity that would be 
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calculated if there were no continuum feature present. 

Many of these emission features are resonance lines 
(R. J. Reynolds 2005, private communication) and, be- 
cause of the substantial IS radiation field at these wave- 
lengths, further work will be required to determine 
whether these lines arise due to collisional excitation in 
the HIM or due to resonant scattering of the IS radiation 
field. We believe that the short wavelength lines are less 
likely to be resonant scatter because the relative inten- 
sity of the s cattered stellar cont inuum in the S band is 
much lower l)Korpela et al.lll998j) . 

Our Civ AA 1548,1551 intensity of 5820±280 LU is 
marginally higher than values for the diffuse ISM re- 
ported by iMartin fc Bowverl 1^990). However Martin 
and Bowyer assumed a featureless continuum rather than 
a continuum with stellar features. If we make the same 
assumption our value drops to 4330±210 LU which is 
more within the range of values they report. 

Our Ovi AA 1032,1038 intensity of 5724±1100 LU is 
somew hat higher tha n the typical value measured by 
FUSE l|Sheltonl 12003]) . This could be a systematic ef- 
fect due to the uncertain calibration of the S band. In 
fact the ratio of the H2 feature intensities between L and 
S suggests that we have underestimated the sensitivity 
of S. Since this is an average over a very large region, it 
is also possible that the Ovi emission is due to several 
small regions with very high intensity, such as SNR. 

The double peaked feature near 1400 A is a blend of 
Siiv AA 1393.8,1402.8 and Oiv] AA 1400.7,1407.4. The 
Oiv] emission is of great interest because, as a semi- 
forbidden line, it is not strongly affected by the radi- 
ation field, but is a direct indicator of collisional pro- 
cesses. Assuming a 2:1 ratio for the Siiv doublet compo- 
nents, our model determines the Oiv] doublet intensity 
to be 1980±300 LU and the Siiv doublet intensity to be 
1430±160 LU. 

We detect high stage ions of two noble gases in our 
spectra, Hen AA 1640 and Nevi AA 999,1005. Although 
the features corresponding to the Nevi lines are statisti- 
cally significant, the poor fit to the underlying H2 bands 
prevents us from placing more than an upper limit to 
these features. Because of the lack of detection of any Oi 
A 1356, we do not expect any contamination of the Hen 
line with Oi A 1641. The Hen line is much brighter than 
would be expected in a solar abundance CIE model (in 
ratio to Cm and Om). This could be due to depletion of 
C and O to ^0.1 solar. A depletion of this magnitude 
would be significantly higher than the accepted values 
for the ISM. 

In Fig. |21 we show equivalent solar abundance emission 
measure (EM) and limits for several species. The results 
are not entirely consistent between species. For exam- 
ple higher EM is necessary to explain the Siiv and Hen 
emission than are compatible with the amount of Cm 
emission. The questionable Nevi detections would indi- 
cate EM 10 x that suggested by the Ovi measurements. 
The bulk of the inconsistency in Nevi is likely to be due 
to undersubtraction of the underlying H2. Much of the 
inconsistency in other lines is likely due to non-CIE ef- 
fects which are present in gas being heated or cooling 



at these temperatures. Some portion of the difference 
could be due to abundance variations (with ISM deple- 
tions of O and C higher than those of Ne, He, and Si), 
sugg 




Tempuroture (K) 

Fig. 3. — Solar abundance emission measure (EM) vs T for 
detected emission lines and significant upper limits. The lines rep- 
resent the loci of EM necessary to produce the observed intensities. 



the accepted ISM values. In future work we will derive 
constraints to non-CIE cooling models and abundances 
using these line ratios. 

5. CONCLUSIONS 

We have presented the FUV spectra of a 30° region 
around the NEP. We detect a variety of emission lines 
from high and low ionization gas. The high ionization 
lines, taken on a per species basis, are consistent with 
EM from 0.001 to 0.005 c m - 6 pc throughout T = 10 4 5 to 
10 5 ' 5 K range. For those few species which have multiple 
lines in the band, the calculated temperatures tend to 
fall near the CIE peak abundance for the species. Despite 
this, a linear combination of CIE models with a restricted 
abundance vs T gradient does not provide a good fit 
to the spectrum. This is likely due to non-CIE effects 
and abundance variations are not directly T related. It 
is also likely that many lines include some resonantly 
scattered stellar radiation. Modeling this scattering will 
be a significant portion of our continued research. 

We have discerned line emission from species that are 
expected to exist in the WNM and WIM. We interpret 
these lines as resonantly scattered stellar radiation. In 
the future, we will investigate if spatial variation of the 
Alii and Sin features can indicate the effect of grain for- 
mation and destruction on the gas phase abundances of 
these elements. 

Despite relatively low galactic N{HI) in this direction, 
we see a substantial amount of H2 fluorescence. Further 
modeling of the H2 fluorescence is necessary to fully 
understand faint spectral features which could be related 
to H2 or due to unrelated atomic lines. 



SPEAR/FIMS is a joint project of KASSI & KAIST 
(Korea) and U.C., Berkeley (USA), funded by the Korea 
MOST and NASA Grant NAG5-5355. 



REFERENCES 



Draine, B. T. 2003, ApJ, 598, 1017 
Edclstein, J., et al. 2005a, ApJ, this volume 



Edelstein, J., et al. 2005b, ApJ, this volume 

Elvis, M., Lockman, F. J., & Fassnacht, C. 1994, ApJS, 95, 413 



SPEAR observations of the NEP 



5 



Knauth, D. C, et al. 2004 BAAS, 205.5906 

Korpela, E. J., Bowyer, S., & Edelstein, J. 1998, ApJ, 495, 317 

Korpela, E. J., et al. 2003, Proc. SPIE, 4854, 665 

Labov, S. E., Yan, B. H. W., & Wai, Y. C. 1989, BAAS, 21, 1123 

Martin, C, & Bowyer, S. 1990, ApJ, 350, 242 

Reed, B. C. 2001, PASP, 113, 537 

Rhee, J. G., et al. 2002, JASS, 19, 57 



Sasseen, T. P., et al. 2002, ApJ, 566, 267 

Savage, B. D., k, Mathis, J. S. 1979, ARA&A, 17, 73 

Shclton, R. L. 2003, ApJ, 589, 261 

Shull, J. M., et al. 2000, ApJ, 538, L73 

Young, P. R., et al. 2003, ApJS, 144, 135 



